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Abstract 

A three-dimensional continuum dislocation theory for single crystals con¬ 
taining curved dislocations is proposed. A set of governing equations and 
boundary conditions is derived for the true placement, plastic slips, and loop 
functions in equilibrium that minimize the free energy of crystal among all 
admissible functions, provided the resistance to the dislocation motion is 
negligible. For the non-vanishing resistance to dislocation motion the gov¬ 
erning equations are derived from the variational equation that includes the 
dissipation function. A simplihed theory for small strains is also provided. 
An asymptotic solution is found for the two-dimensional problem of a single 
crystal beam deforming in single slip and simple shear. 

Keywords: dislocations (A), crystal plasticity (B), hnite strain (B), 
variational calculus (C) . 


1. Introduction 


In view of a huge number of dislocations appearing in plastically deformed 
crystals (which typically lies in the range 10® 10^® dislocations per square 

meter) the necessity of developing a physically meaningful continuum dis¬ 
location theory (CDT) to describe the evolution of dislocation network and 
predict the formation of microstructure in terms of mechanical and thermal 
loading conditions becomes clear to all researchers in crystal plasticity. One 
of the main guiding principles in seeking such a continuum dislocation the¬ 
ory has hrst been proposed by Hansen and Kuhlmann-Wilsdorf (|1986 ) 
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form of the so-called LEDS-hypothesis: the true dislocation structure in the 
hnal state of deformation minimizes the energy of crystal among all admis¬ 
sible dislocation conhgurations. In view of numerous experimental evidences 
supporting this hypothesis (see, i.e., (Hughes and Hansen, 1997 Kuhlmann- 


Wilsdorf, 1989, 2001 Laird et al., 1986)), its use in constructing the con¬ 
tinuum dislocation theory seems to be quite reasonable and appealing. For 
the practical realization one needs to i) specify the whole set of unknown 
functions and state variables of the continuum dislocation theory, and ii) lay 
down the free energy of crystals as their functional to be minimized. Such 
program has been implemented by Berdichevsky (2006a) in the linear, and by 


Le and Gunther (2014) in the nonlinear setting of CDT for networks of dislo¬ 


cations, whose lines are straight and remain so during the whole deformation 
process (see also (Ortiz and Repetto, 1999; Ortiz et ah, 2000)). The devel¬ 


oped CDT has been successfully applied to various two-dimensional problems 
of dislocation pileups, bending, torsion, as well as formation of dislocation 


patterns in single crystals (see 1 

Berdichevsky and Le 

2007; 

Kaluza and Le 

20TT| 

Kochmann and Lei |2008a 


2009[ 

Foster et al. 

2015; 

Le and Sembir- 


ing, 2008a|[b , 2009[ Le and Nguyenf 2012 2013)). Let us mention here the 


similar approaches suggested in (|Acharya and Bassani 2000[ Acharya, 2001| 


Engels et al.l |20121 |GurtinH2002t [Gurtin et al.[|2007[ [Mayeur and McDowell 


2014 


Oztop et ah, 2013) which do not use the LEDS-hypothesis explicitly 


but employ instead the extended principle of virtual work for the gradient 
plasticity. However, as experiences and experiments show, dislocation lines 
are in general loops that, as a rule, can change their directions and curva¬ 
tures depending on the loading condition and crystal’s geometry. Therefore 
the extension of CDT to networks of dislocations whose lines are curves in 
the slip planes is inevitable. To the best of author’s knowledge, such three- 
dimensional continuum dislocation theory based on the LEDS-hypothesis for 
curved dislocations has not been developed until now. It became also clear 
to him that the latter’s absence was due to the missing scalar dislocation 
densities for the network of curved dislocations. 

The hrst attempt at constructing a continuum theory that can predict in 
principle not only the dislocation densities but also the direction and curva¬ 
ture of the dislocation lines has been made by Hochrainer et al. (2007) in form 
of the so-called continuum dislocation dynamics. Their theory starts with 
the dehnition of the dislocation density that contains also the information 
about the orientation and curvature of the dislocation lines. Then the set of 
kinematic equations is derived for the dislocation density and curvature that 
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requires the knowledge about the dislocation velocity. The relation between 
the dislocation density and the macroscopic plastic slip rate via the disloca¬ 
tion velocity is postulated in form of Orowan’s equation. The couple system 
of crystal plasticity and continuum dislocation dynamics becomes closed by 
the constitutive equation of a flow rule type (see (Hochrainer et ah, 2014 


Sandfeld et ah, 2011, 2015 Wulhnghoff and Bohlke, 2015[ )). In addition to 


the heavy computational cost of such theory, the relation to thermodynam¬ 
ics of crystal plasticity and to the LEDS-hypothesis is completely lost: the 
equilibrium solution found in this theory may not minimize the energy of 
crystal among all admissible dislocation conhgurations. Let us mention also 


a continuum approach proposed recently by Zhu et ah (2013); Zhu and Xiang 


(2014) in which the three-dimensional dislocation structure is characterized 
by two families of disregistry functions that may take only integer values. The 
dislocation density can then be expressed in their terms. The coupled system 
of equations is derived from the underlying discrete dislocation dynamics for 
the displacement and disregistry functions. This approach is subject to the 
same critics as that proposed in (Hochrainer et ah, 2007). 

The aim of this paper is to extend the nonlinear continuum dislocation 
theory (CDT) developed recently by Le and Giinther (2014) to the case of 
crystals containing curved dislocations. Provided the dislocation network is 
regular in the sense that nearby dislocations have nearly the same direction 
and orientation, we introduce a loop function whose level curves coincide 
with the dislocation lines. Taking an inhnitesimal area perpendicular to the 
dislocation line at some point of the crystal, we express the densities of edge 
and screw dislocations at that point through the resultant Burgers vectors of 
dislocations whose lines cross this area at right angle. Such scalar densities 
contain not only the information about the number of dislocations, but also 
the information about the orientation and curvature of the dislocation lines. 
In case of dislocation motion we introduce the vector of normal velocity of 
dislocation line through the time derivative of the loop function. Following 


Kroner (1992) and (Berdichevsky, 2006b) we require that the free energy 


density of crystal depends only on the elastic strain tensor and on the above 
scalar densities of dislocations. Then we formulate a new variational princi¬ 
ple of CDT according to which the placement, the plastic slip, and the loop 
function in the hnal state of equilibrium minimize the free energy functional 
among all admissible functions. We derive from this variational principle 
a new set of equilibrium equations, boundary conditions, and constitutive 
equations for these unknown functions. In case the resistance to dislocation 
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motion is significant, the variational principle must be replaced by the vari¬ 
ational equation that takes the dissipation into account. The constructed 
theory is generalized for single crystals having a hnite number of active slip 
systems. We provide also the simplihcations of the theory for small strains. 


As compared to the continuum dislocation dynamics proposed in (Hochrainer 


et ah, 2007 Zhu and Xiang, 2014) our theory is advantageous not only in the 


computational cost due to its simplicity, but also in its full consistency with 
the LEDS-hypothesis. In the problem of single crystal beam having only one 
active slip system and deforming in simple shear, the energy minimization 
problem reduces to the two-dimensional variational problem. We solve this 
problem analytically for the circular cross section and asymptotically for the 
rectangular cross section. We will show that this solution reduces to that 


found in (Berdichevsky and Le, 2007) for the crystals with thin and long 


cross-section. 

The paper is organized as follows. After this short introduction we present 
in Section 2 the three-dimensional kinematics for single crystals deforming 
in single slip. Section 3 formulates the variational principles of the three- 
dimensional CDT and derives its governing equations. Section 4 extends 
this nonlinear theory to the case of single crystals with n active slip sys¬ 
tem. Section 5 studies the three-dimensional small strain CDT. Section 6 is 
devoted to the analytical and asymptotic solutions of the two-dimensional 
energy minimization problem of a single crystal beam deforming in simple 
shear. Finally, Section 7 concludes the paper. 


2. 3-D kinematics for single crystals deforming in single slip 

Nonlinear CDT starts from the basic kinematic resolution of the defor¬ 


mation gradient F = dy/dx into elastic and plastic parts (Bilby et al., 1957) 

F = • FP. (1) 

We attribute an active role to the plastic deformation: F^ is the deformation 
creating dislocations (either inside or at the boundary of the volume element) 
or changing their positions in the crystal without distorting the lattice par¬ 
allelism (see Fig. [^. On the contrary, the elastic deformation F® deforms 
the crystal lattice having frozen dislocations (Le and Gunther, 2014). Note 


that the lattice vectors remain unchanged when the plastic deformation is 
applied, while they change together with the shape vectors by the elastic 
deformation. 
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Figure 1: Multiplicative decomposition 

We consider first a single crystal deforming in single slip. In this case 
let ns denote the right-handed triad of nnit lattice vectors of the active slip 
system by s, p, and m, where s points to the slip direction, p lies in the slip 
plane and is perpendicnlar to s, and m is normal to the slip plane. Withont 
restricting generality we may choose the rectangnlar cartesian coordinate 
system {xi,X 2 ,xz) in the reference confignration snch that its basis vectors 
coincide with these lattice vectors (see Fig. 

ei = s, e 2 = p, 63 = m. 

The plastic deformation is then given by 

= I +/3(x)s (g) m = I-f/3(x)ei (g) es, (2) 

with 13 being the plastic slip. We assnme that all dislocations cansing this 
plastic deformation lie completely in the slip planes and the dislocation net¬ 
work is regnlar in the sense that nearby dislocations have nearly the same 
direction and orientation. This enables one to introdnce a scalar fnnction 
l{xi,X 2 ,X 3 ) (called a loop fnnction) snch that its level cnrves 


/(a;i,a; 2 ,C 3 ) = c, (3) 

with C 3 and c being constants, coincide with the dislocation lines. Thns, 
in this three-dimensional kinematics we admit, according to eqnation (|^, 
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only the conservative motion of dislocations and exclude from considera¬ 
tion the dislocation climb which is an important mechanism of temperature- 
dependent creep. We denote by v and r the plane unit vectors normal and 
tangential to the dislocation line. From equation ([^ follow 

^ = —I (^,1^1 + ^,2^2), T = — { — 1 , 2^1 + ^,1^2), 

vh+F \/F+F 

where the comma before an index denotes the partial derivative with respect 
to the corresponding coordinate. Note that r, m form a right-handed 
basis vectors of the three-dimensional space (see Fig. |^. 



Figure 2: A dislocation loop in the chosen coordinate system 


Ortiz and Repetto (1999) introduced the resultant Burgers vector of ex¬ 


cess dislocations, whose lines cross the area A in the reference conhguration, 
in the following way 


br = ® • dx, 


( 4 ) 


where C is the close contour surrounding A. Le and Giinther (2014) have 


shown that, in the continuum limit, when the atomic distance goes to zero 
at the fixed sizes of the representative volume element and the fixed density 
of dislocations per area of unit cell, integral (|^ gives the total closure failure 
induced by which must be equal to the resultant Burgers vector. It is 
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natural to assume continuously differentiable in this continuum limit, so, 
applying Stoke’s theorem we get from Q 

b, = - / (FP X V) ■ nda, 

Ja 

where x denotes the vector product, da the surface element, and n the unit 
vector normal to A. This legitimates the introduction of the dislocation 
density tensor 

T = -FP X V. 

For the plastic deformation taken from ([^ 

T = -F^' X V = s (g) (V/5 X m). 

If we choose now an inhnitesimal area da with the unit normal vector r, then 
the resultant Burgers vector of all excess dislocations, whose dislocation lines 
cross this area at right angle is given by 

br = T ■ r da = —s (V/d ■ iy)da = —s 8^(5 da. 

This resultant Burgers vector can be decomposed into the sum of two vectors 


br = br± + b^ii = —{usy + TSr)dyf3 da, 


where = s ■ u = ui and = s ■ r = ri are the projections of the 
slip vector onto the normal and tangential direction to the dislocation line, 
respectively. This allows us to dehne two scalar densities (or the numbers of 
excess dislocations per unit area) of edge and screw dislocations 


P± = 


P\\ = 



^\sudu/3\ 

^Isrdu/Sl 


1 

b 

1 

b 


+ 1^,2^,2) 


(5) 


with b the magnitude of Burgers vector. We see that the three-dimensional 
dislocation densities p± and p\\ depend on both the gradient of the plastic 
slip and the gradient of the loop function Z(x) through the vectors u and r. 

Consider now the case of motion of dislocation loops in the slip plane. In 
this case we allow the loop function to depend explicitly on time t such that 
equation 

l{xi,X 2 ,C 3 ,t) = C (6) 
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with fixed constants C 3 and c describes one and the same dislocation line 
during its motion in the slip plane. Letting ( be the variable along the dislo¬ 
cation line, we may represent the level curve defined by ([^ in the parametric 
form 

Xi = Xi(C,t), X2=X2iC,t). 

Fixing ( and taking the differential of (|^ we obtain 

l ^diXi l 2dx2 — 0 


that yields 


dxi dx2 


Since v = ^ei -|- ^^2 is the velocity of the fixed point on the dislocation 
line with coordinate (, we define the normal velocity of the dislocation line 
as follows 


Vu = W -U = 



dxi 
^ dt 



i 



(7) 


with I = This kinematic quantity will be used in the case of non-vanishing 
resistance to dislocation motion. 


3. Governing eqnations for single crystals deforming in single slip 


According to Kroner (1992), the elastic deformation F® and the dislo¬ 


cation densities p± and py characterize the current state of the crystal, so 
these quantities are the state variables of the continuum dislocation theory. 
The reason why the plastic deformation cannot be qualified for the state 
variable is that it depends on the cut surfaces and consequently on the whole 
history of creating dislocations. Likewise, the gradient of plastic strain ten¬ 
sor cannot be used as the state variable by the same reason. In contrary, 
the dislocation densities depend only on the characteristics of dislocations 
in the current state (Burgers vector and positions of dislocation lines) and 
not on how they are created, so p± and py, in addition to F®, are the proper 
state variables. Thus, if we consider isothermal processes of deformation, 
then the free energy per unit volume of crystal (assumed as macroscopically 
homogeneous) must be a function of F®, p±, and py 


V’ = V'(F^p±,py). 
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Now, if we superimpose an elastic rotation R onto the actual deformation of 
the body, then the total and elastic deformation change according to 

F* = R F, F^* = R Fh 


At the same time, the dislocation densities p_[_ and p\\ remain unchanged. 
As such superimposed elastic rotation does not change the elastic strain and 
the dislocation densities, we expect that the energy remains unchanged. The 
standard argument (see, e.g., (Gurtin, |l981 )) leads then to 

V' = '0(C",P±,P||), 


where C® is the elastic strain dehned by 

C® = F*^^ ■ F® 


Let the undeformed single crystal occupy some region V of the three- 
dimensional space. The boundary of this region, dV, is assumed to be the 
closure of union of two non-intersecting surfaces, dk and dg- Let the place¬ 
ment be a given smooth function of coordinates 

y(x) = X uo(x) at dk, (8) 

where Uo(x) = y(x) — x is the given displacement vector. Such condition 
does not admit dislocations to reach this part dk of the boundary, so we set 

/3(x) = 0, /(x) = 0 at dk- (9) 

At the remaining part dg the “dead” load (traction) t is specihed. Note that, 
in case the whole boundary is free, we do not have any kinematic constraint 
at dV. However, as the energy density is invariant with respect to the shift 
of the loop function /(x) on an arbitrary constant which does not change 
the dislocation densities, we can impose on this scalar function the following 
constraint 

/ /(x) dx = 0, 

Jv 

where dx = dxidx 2 dx^ denotes the volume element. If no body force acts on 
this crystal, then its energy functional is dehned as 

/[y(x),/3(x),/(x)] = f w{F,(3,V(3,Vl)dx - f t ■ y da, (10) 
Jv Jds 
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where 


w;(F,AV/3,V0 = ^(C^px,Pll), 


( 11 ) 


Provided the resistance to the dislocation motion is negligibly small and no 
surfaces of discontinuity occur inside crystals, then the following variational 
principle is valid for single crystals with one active slip system: the true place¬ 
ment vector y(x), the true plastic slip /3(x), and the true loop function /(x) 


in the final equilibrium state of deformation minimize energy functional (10) 


among all continuously differentiable helds y(x), /3(x), and /(x) satisfying 
constraints ([^ and (|^. 

Let us derive the equilibrium equations from this variational principle. 


We compute the hrst variation of functional (10) 


61 = 


dw dw dw \ f 


where P = dw/dF. Integrating the hrst, third, and fourth term by parts 
with the help of Gauss’ theorem and taking the conditions ([^ and ([^ into 
account, we obtain 


61 = / [-5y ■iP-V) + {w 0 -V- w^p) 6 fi - (V ■ w^i) 6 l] dx 


+ / [(P ■ n — t) ■ 5y -|- tcv /3 ■ n 5/3] da -|- / Wyi ■ n 6 l da = 0. (12) 
Jds JdV 


Equation (12) implies that the minimizer must satisfy in V the equilibrium 
equations 

P ■ V = 0, -Wi3 + V ■ tcv/? = 0, V ■ tcv/ = 0, (13) 

subjected to the kinematic boundary conditions (|^ and (§ at dk, and the 
following natural boundary conditions 

P • n = t, tav /3 ■ n = 0, w\/i ■ n = 0 at c/s- (14) 

We call P the hrst Piola-Kirchhoh stress tensor, = —w^ the resolved 
shear stress (or Schmid stress), and ? = —V ■ wv /3 the back stress. The hrst 


equation of (13) is nothing else but the equilibrium of macro-forces acting on 


the crystal, the second equation represents the equilibrium of micro-forces 
acting on dislocations, while the last one expresses the equilibrium condition 
for the curved dislocation lines. 
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The constitutive equations for P = w-p, Wvi 3 i and w^i can easily be 
obtained from the free energy density (0. First, we express F® in terms of 
F and (3 with the use of Q and ([^ 

F" = F • FP-^ = F ■ (I - /3s ® m). 


Now, the standard differentiation using the chain rule and the above relation 
yields the hrst Piola-Kirchhoff stress tensor 


P = m;p = 2F" ■ V’c- ■ 

For the resolved shear stress (Schmid stress) we get 

Tr = —Wjj = 2s ■ F^“^ ■ C® • ■ m. 


(15) 


(16) 


Likewise, from (|^ follows 

^v/3 = ^ sign(sj.a^/3)s^ + ?/>p||Sign(s^a^/3)s^]j/. (17) 

Thus, the vector wy/? is two-dimensional. Finally, we compute Wyz directly 
in components using formulas (|^. Since p±_ and p\\ do not depend on 
so wi^ = 0, and the vector w^/i is also two-dimensional. For its hrst two 
components we have 


WZi = 


tl^p^siga^Sud^P) 

+ ^jJp^^sign{sr^^(3) 
wi^ = ]- 'ipp^siga{s^d^/3) 


2/2(/3,i/,i + /3,2/,2) 2/3,i/,i-7/3,2/,2 


+ i%y i‘^i + 1% 

2 ^, 1 ^, 2 ( Ai ^,1 + A 2 ^, 2 ) Ai ^,2 


ill + 1 %? 


II + II 


2/ J 


‘^lpipi!^pl ,1 + f^plp) t^plp 


ill+ll? 


ll + 1 % 


-7 i)p sign{srdy(3) 


2/2(/3,iA+A2/,2) AA,1 + 2A2^ 


{ll+Ilf 


+ 


Ip + Ip 


(18) 


Substituting the constitutive equations (15)-(18) into (13)-(14) we get the 
completely new system of equations and boundary conditions which, together 
with ([^ and 0. enable one to determine y(x), /3(x), and /(x). Note that 
equations (13)i and (13)2 are coupled via the hrst Piola-Kirchhoh stress 
tensor and the Schmid stress containing both F and (3, while equations (13)2 
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and (13)3 are coupled because both contain the gradients of /3 and 1. All 


equations are strongly nonlinear partial differential equations. 

The above theory has been developed for the case of negligibly small resis¬ 
tance to dislocation motion and plastic slip. In real crystals there is however 
always the resistance to the dislocation motion and plastic slip causing the 
energy dissipation that changes the above variational principle as well as the 
equilibrium conditions. We assume that the dissipation function depends on 
the plastic slip rate f3 and on the normal velocity of the dislocation loop 
given by ([^ (or, equivalently, on /). Thus, 

D = D0J), 

When the dissipation is taken into account, the above formulated variational 


principle must be modihed. Following (Sedov, 1965 Berdichevsky, 1967) we 


require that the true placement y(x, t), the true plastic slips /3(x, t), and the 
true loop function /(x, f) obey the variational equation 


61 + (— -6/3 H- -61) dx = 0 

Jv d/3 dl 


(19) 


for all variations of admissible helds y(x, f), /3(x, f), and /(x, t) satisfying the 
constraints ([sj) and §. Together with the above formula for 61 and the 


arbitrariness of 6 y, 6/3, and 61 in V as well as at dg, equation (19) yields 


P-V = 0, 


_ dD 

-W,3 + V-Wyp = 

d/3 


V • ws/i - (20) 

dl 


which are subjected to the kinematic boundary conditions ([^ and ([^, and 
the natural boundary conditions (14). The constitutive equations remain ex¬ 


actly the same as (15)-(18). For the rate-independent theory the dissipation 


function can be assumed in a simple form 

D = K,\^\ + K2\i\, 

with Ki and K 2 being positive constants. We call Ki the critical resolved 


shear stress and K 2 the Peierls threshold. In this case equations ( 20 ) 2,3 
become 

-Wj 3 + V ■ tcv/? = Ki sign/3, V ■ w^/i = K 2 sign/ 

for non-vanishing /3 and /. These are the yield conditions for /3 and /: (3 and 
I are non-zero if and only if 


W/3 + V • Wv/?! = A'l, \V-w^i\=K2 . 
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On the contrary, if the expressions on the left-hand sides are less than Ki 
and i^ 2 , the plastic slip cannot evolve and the dislocation lines cannot move: 
/3 = 0 and I = 0. Thus, they are frozen in the crystal. 


4. Extension to multiple slip 

The extension to the case of single crystals having n active slip systems 
can be done straightforwardly under the assumptiorj^ 

n 

= I + ^/3“(x)s“®m“, (21) 

a=l 

with /3“ being the plastic slip, where the pair of constant and mutually or¬ 
thogonal unit vectors s“ and m“ is used to denote the slip direction and the 
normal to the slip planes of the corresponding a-th slip system, respectively. 
Here and later, the Gothic upper index a running from 1 to n numerates the 
slip systems, so one could clearly distinguish /3“ from the power function. We 
denote by p“ the unit vector lying in the slip plane such that s“, p“, and m“ 
form a right-handed basis vectors. For each slip system we can introduce the 
coordinates associated with these basis vectors 


e“ = s“-x, e2 = p“-x, e3“ = m“-x. 


( 22 ) 


Equations (22) can be regarded as the one-to-one linear transformation re¬ 


lating and X according to 

= M“ ■ X, X = 

where M“ is the 3x3 matrix whose rows are basis vectors s“, p“, and m“. 
Thus, any function of x can be expressed as function of and vice versa. 
For the plastic slip caused by dislocations of the slip system a we assume 
that their lines lie completely in the slip planes parallel to the (^“, ^ 2 )"Plcme. 
To describe the latter we introduce the loop function '^ 2 )■Cs) such that 
its level curves 

nei,e2,c3) = c, ( 23 ) 


^In conventional crystal plasticity the kinematic equation for is usually formulated 


in rate form that does not always reduces to (211 (Ortiz and Repetto, 1999 Ortiz et al. 


2000 ). 
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where C 3 and c are constants, coincide with the dislocation lines. We denote 
by i/" and t“ the plane unit vectors normal and tangential to the dislocation 
line. From equation (23) follow 

1 


= 


(^“ 1 )^ + (1%) 


:(/“s“ + r 2 p“), 


t“ = 


1 


(KiY + (1%) 




where the semicolon in indices denotes the partial derivatives of the loop 
function with respect to so these vectors lie in the slip planes parallel 

to the (^ 1 , ^ 2 )-plane as expected. 

For the plastic deformation (21) the dislocation density tensor becomes 

n 

T = -FP X V = ^ s“ O (V/3“ X m“). (24) 

a=l 

To characterize the geometrically necessary dislocations belonging to one slip 
system we consider one term T“ = s“ ( 8 ) (V/3° x m“) in the sum (24). Let us 
choose an inhnitesimal area da with the unit normal vector t“ and compute 
the resultant Burgers vector of all excess dislocations of the system a, whose 
dislocation lines cross this area at right angle 

b“ = T“ ■ t“ da = —s“ (V/9“ ■ iy'^)da = —s“ d^a(3° da. 

This resultant Burgers vector can be decomposed into the sum of two vectors 

b“ = b“_L + b“|| = + T‘'s°)d^a(3^ da, 

where = s“ ■ and s“ = s" • are the projections of the slip vector onto 
the normal and tangential direction to the dislocation line, respectively. This 
allows us to dehne two scalar densities of edge and screw dislocations of the 
corresponding slip system 


p 1 = 


lb“ I 
b 

|b“ I 


P\\ 


= = I 




imiKi+ 


+ (/“ )2 

^;2(AV“i+AV“2) 


+ (^“ 2 ) 


(25) 


We see that the dislocation densities and py depend only on the partial 
derivatives (3°’^ and a = 1, 2. For the moving dislocations we allow the 
loop functions to depend on time t such that the level curves 


r{ei,C2,C3,t)=C, 
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with C 3 and c being constants, coincide with the dislocation lines dnring 
their motion. Similar to the single slip we introdnce the normal velocities of 
dislocation lines as follows 

/“ 

with I = Zt. These kinematic qnantities will be nsed in the model with 
dissipation. 

From the above discnssion of kinematics we see that a single crystal with 
n active slip systems is a generalized continnnm with 3 + 2n degrees of free¬ 
dom at each point: y(x), /?“(x), and Z“(^“(x)), a = l,...,n. We reqnire 
that the free energy per nnit volnme of crystal (assnmed as macroscopically 
homogeneons) mnst be a fnnction of C® = ■ F® (where F® = F ■ F^“^), 

and p]| 

Under the same loading condition as for the crystal with single slip we write 
down the energy fnnctional 



J[y(x),/3“(x),Z“(r(x))]= / / t-yda, (26) 

Jv Jds 

where 

n;(F,/?^V/?^V/“) = ^(C^pl,p^). 

Provided the resistance to the dislocation motion is negligibly small, we for¬ 
mulate the following variational principle for single crystals with n active slip 
systems: the true placement vector y(x), the true plastic slips /3“(x), and 
the true loop functions /“(^“(x)) in the final equilibrium state of deformation 
minimize energy functional (26) among all continuously differentiable helds 
y(x), /3“(x), and Z“(^“(x)) satisfying the constraints 

y(x) = X -F uo(x), /3“(x) = 0, Z“(x) = 0 at dk, (27) 


Applying the same calculus of variation and taking into account the ar¬ 
bitrariness of the variations of y(x), /3(x), and Z(^“(x)) in V as well as at dg, 
one can show that the minimizer must satisfy in V the equilibrium equations 


P • V = 0, —ui^ + V ■ rcv/ 3 “ = 0, V • = 0, (28) 
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subjected to the kinematic boundary conditions (27) at dk and the following 
natural boundary conditions 


P ■ n = t, tcv/ 3 “ ■ n = 0, ■ n = 0 at (29) 

The constitutive equations for P = tap, —Wjsa, w^/pa, and w^ia can easily be 
obtained from the above free energy density by standard differentiation. For 
the hrst Piola-Kirchhoff stress tensor and the Schmid stresses we have 


p = tof = 2F' ■ V'c- ■ 
r; = -wp. = 2s" • F'’-'^ • C“ ■ V'c- ■ m". 


(30) 

(31) 


Likewise, from (25) follows 




(32) 


Thus, the vectors wv/ 3 “ are two-dimensional. Finally, for we have 


/ dp‘', op\\\ ( op\\\ 

»v,. = (<^,1 ^ S" + (^Vl ^ P“. (33) 

Differentiating formulas (25) for the dislocation densities p']_ and py with 
respect to Z" and l° 2 , we get 


1 

^ = ^sign(4a.«/3“) 

dpn 1 

— = -sign(s“a,«/3“) 
^ = ^sign(s;9„.;3“) 

dpu 1 

— = -sign(4a,«/3“) 


2(/“i)nAV.u+Ay“2) 2/3 “Z“i + AV“2 

+ (^“ 2 )^ 


2^;1^;2(AV“i+AV“2) 


2'\2 


+ 


(itir + (^“ 2 )^ 


+ 




my + m. 

2/;l/MAV;l + 

\iim+ii%yy ' im+iim 
aim + iim+nm 


Substituting the constitutive equations (30)-(33) into (28)-(29) we get the 
completely new system of equations and boundary conditions which, together 
with (27), enables one to determine y(x), /3“(x), and Z“(x). 
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For the case of non-zero resistance to dislocation motion leading to the 
energy dissipation we take the dissipation function in the form 

D = D0\n, 

We require that the true placement y(x, t), the true plastic slips /3“(x, t), and 
the true loop functions Z“(x, t) obey the variational equation 


51 + 


[ f2{^5(3^ + dx = 0 


(34) 


for all variations of admissible helds y(x, t), /3“(x, f), and /“(x, t) satisfying 
the constraints (27). It is then easy to show by exactly the same arguments 
like those used at the end of the previous Section that equation (34) yields 

dD 


P ■ V = 0, —Wpa + V ■ tCV/3“ = 




_ dD 

V ■ Wv/o — 7^, 


which are subjected to the boundary conditions (27) and (29). The consti¬ 
tutive equations remain exactly the same as (30)-(33). 


5. Small strain theory 


Let us simplify the above theory for small strains. In this case, instead of 
the placement y(x) we regards the displacement u(x) that is related to the 
former by 

u(x) = y(x) - X 

as the unknown function. Thus, the total compatible deformation is 


F 


dx. 


I + uV. 


We assume that the displacement gradient uV (called distortion) is small 
compared with I. Concerning the plastic deformation given by (21) we also 
assume that the plastic slips (5° are much smaller than 1. Using the the 
multiplicative resolution ([^ to express F® through F and F^“^ and neglecting 
the small nonlinear terms in it, we get 


F^ = F ■ F^-^ = I + /3" 
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where 

n 

/3® = uV — /3“s“ ® m“. 

a=l 

The last equation can be interpreted as the additive resolution of the total 
distortion into the plastic and elastic parts. The total compatible strain 
tensor held can be obtained from the displacement held according to 

e = ^(uV + Vu). 

The incompatible plastic strain tensor held is the symmetric part of the 
plastic distortion held 

1 1 

= -(/3 + 0^) = - (g) m“ + m“ ® s“). 

a=l 

Accordingly, the elastic strain tensor held is equal to 


The dislocation densities remain exactly the same as in the hnite strain 


theory. They are given by the formulas (25). Concerning the free energy 


density we will assume that it depends on the elastic strain e® and on the 
dislocation densities p\ and py 

The energy of crystal containing dislocations reads 


/[u(x),/3(x),/(x)] = / t(;(uV,/3“, V/9“, V/“) dx — / t-uda, (35) 
Jv Jds 

where tc(uV,/?“, V/3“, V/“) = p±, p\\). Provided the resistance to the 

dislocation motion can be neglected, then the following variational principle 
is valid for single crystals: the true displacement held u(x), the true plastic 
slips /5“(x), and the true loop functions P(x) in the final state of deformation 


in equilibrium minimize energy functional (35) among all admissible helds 
satisfying the constraints 


u(x) = uo(x), /3“(x) = 0, /“(x) = 0 at dk- 


(36) 
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The standard calculus of variation similar to the previous case leads to 
the equilibrium equations 


cr • V = 0, tCfl — V ■ tCv/ 3 “ =0, V ■ w^ia = 0, 


(37) 


subjected to the kinematic boundary conditions (36) at dk and the following 
natural boundary conditions 

cr • n = t, tCv/ 3 “ ■ n = 0, w^ia ■ n = 0 at 5^. (38) 

The constitutive equation for the Cauchy stress tensor becomes 

d%l) 


(j = 


de^ 


(39) 


For the Schmid stress we obtain 

r“ = —wpa = s" • cr • m“. (40) 

The constitutive equations for tcv/ 3 “ and w^ia remain unchanged as com¬ 
pared with (32) and (33). The new set of governing equations and boundary 
conditions are obtained by substituting (39), (40), (32), and (33) into the 
equilibrium equations (37) and boundary conditions (38). Note that even for 
small strain the system of governing equations remain as a whole nonlinear. 

It is a simple matter to modify the theory for the case of non-zero resis¬ 
tance to dislocation motion and plastic slip leading to the energy dissipation. 


6. Simple shear deformation of a single crystal beam 

Let us consider now the simple shear deformation of a single crystal beam 
having only one active slip system. The crystal occupies in its initial conhgu- 
ration a long cylinder of an arbitrary cross section such that (xi, X 2 ) € A and 
0 < X 3 < L (see Fig. [^for the beam of rectangular cross section). As before, 
the slip system is chosen such that the vectors s, p, and m coincide with ei, 
62 , and 03 , respectively. We realize the simple shear deformation by placing 
this crystal beam in a “hard” device with the prescribed displacements at 
the boundary of the crystal such that 

yi = xi + 7 x 3 , y 2 = X 2 , ys = X 3 . (41) 

We assume that the length of the crystal L is large enough compared with the 
sizes of the cross section to guarantee the uniform simple shear deformation 
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Figure 3: Simple shear deformation of a beam of rectangular cross section 


state for which equation (41) is valid everywhere in the beam. If the overall 


shear 7 is sufficiently small, then it is natural to expect that the crystal 
deforms elastically and the plastic slip as well as the loop function vanish. If 
this parameter exceeds some critical threshold, then dislocation loops may 
appear (see one dislocation loop in Fig. |^. Due to the almost translational 
invariance in xs-direction we may assume that /3(x) and /(x) depend only on 
xi and X 2 (except perhaps the neighborhoods of X 3 = 0 and X 3 = L). 

For simplicity let us consider the small strain theory. Then the only non¬ 


zero components of the total strain tensor, under the condition that (41) is 
valid everywhere, are 

1 

— ^31 — 27 - 

Since the plastic distortion tensor (3 has only one non-zero component /Sia = 
I 3 {xi,X 2 )-, the non-zero components of the elastic strain tensor read 


£13 = £31 = 1^(1 - 


(42) 


With the loop function being l{xi,X 2 ) the dislocation densities are given by 


P± 

P\\ 


1 

b 

1 

b 


+ / 3 , 2 ^, 2 ) 

+ 13,21,2) 

^7 + ^7 


(43) 
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For the small strain theory we propose the free energy per nnit volnme 
of the nndeformed crystal in the most simple form 




^A(tr£")2 + /itr(£" ■ e^) + 

^ ^ P S ^ Ps 


(44) 


Here e® is the elastic strain tensor, A and are the Lame constants, ki and 
p 2 are material constants, while ps can be interpreted as the satnrated dis¬ 
location density. The hrst two terms in (44) represent the free energy of 
the crystal dne to the macroscopic elastic strain, where we assnme that the 
crystal is elastically isotropic. The last two terms in (44) correspond to the 
energy of the dislocation network for moderate dislocation densities (Gnrtin| 


2002 Gnrtin et ah, 2007). Note that, for the small or extremely large dislo¬ 


cation densities close to the satnrated valne, the logarithmic energy proposed 
by Berdichevsky ( |2006b|[a ) tnrns ont to be more appropriate. Snbstitnting 
formnlas (42) and (43) into (44) we obtain 


«;(/?, V/?, V/) = ^/i(7 






2h^pl L 


ki 




ill + 1%? 

+ /C 2 


+ f^pipY 

ill+ 11? 


Since in this case the side bonndary does not allow dislocations to reach it, 
we can pose on both fnnctions (5{xi,X2) and l{xi,X 2 ) the Dirichlet bonndary 
conditions 


(3{xi,X2) = ^, lixi,X2) = ^ ior {xi, X2) e dA, 


(45) 


The variational problem rednces to minimizing the two-dimensional fnnc- 
tional 

I[/3{xi,X2),lixi,X2)] = L / w{/3,'\//3,'\/l)dxidx2 (46) 

Ja 

among all admissible fnnctions /3(a;i, X 2 ) and /(a;i, X 2 ) satisfying the bonndary 
conditions (45). 

It is convenient to simplify the fnnctional and minimize it in the dimen¬ 
sionless form. Introdncing the dimensionless variables and qnantity 


Xi = hpsXi, X2 = bpsX2, ixi,X2) E A, I = 


pL ’ 
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we can write functional (46) in the form 


I = - 

2 


U L 


(7 - + h 


+ A 2 ^, 2 )^ + (^,21,2) 


m+iiY 


+ k2 


m+iiY 


dxidx2, 
(47) 

where the bar over the quantities are dropped for short. The problem is 
to minimize functional (47) among all admissible functions P{xi,X 2 ) and 
l{xi,X 2 ) satisfying Dirichlet boundary conditions (45). Since the integrand 
of (47) is positive dehnite, the existence of the minimizer in this variational 
problem is guaranteed. 

In one special case the problem degenerates and admits an analytical 
solution. Indeed, if we choose ki = k 2 = k and take ^ to be a circular 
cross section whose boundary is given in the polar coordinates by r = i?, 
then due to the symmetry we may assume that both f3 and I are functions 
of r = ^Jx‘l + X 2 only. It is now a simple matter to show that 

+ A2A)' A 2 (AiA + A2A)' _ .2 

T /70 . 70 \ r, P r' 


(A+A 2 )' 


(A+A 2 )' 


Thus, functional (47) (normalized by 27r) does not depend on l{r) and takes 
the form 

I=f [lil - + lkf3f^]rdr, 


which leads to the Euler equation 


k 


h- ^) + -iYrr),r = 0 . 

This is nothing else but the inhomogeneous modihed Bessel equation that 
yields the following solution (regular at r = 0 ) 

f3[r) = 7 + CIo{r/y/k), 

with Io{x) being the modihed Bessel function of the hrst kind. The coefficient 
C must be found from the boundary condition P{R) = 0 giving 


Thus, 


C = —f/Io{R/Vk). 
/3{r)=^[l-Io{r/Vk)/IoiR/Vk)]. 


( 48 ) 
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The plot of solution (48) for different values of 7 is shown in Fig. where 
we choose the radius such that R = hpgR = 1 and the parameter k = 10“^. 
As 7 increases, the plastic slip increases too. It is seen from this Figure 
that the plastic slip is nearly constant in the middle of the cross section and 
changes strongly only in the thin layer in form of ring near the boundary. 
Since I is a function of r, dislocations in form of circles pile up against the 
boundary of the cross section, leaving the middle of the cross section almost 
dislocation-free. Since the dislocation loops are circles, they have the purely 
edge character at (p = 0 and p = n and the purely screw character at (p = 
± 7 r/ 2 . For all other angles the dislocation loops have the mixed character. 
Note that if ki 7 ^ ^ 2 , fhe strictly axi-symmetric solution is no longer valid 
because the contributions of the edge and screw components to the energy of 
dislocation network are not equal. Note also that, if the logarithmic energy 
is used instead of the quadratic energy, f3 becomes non-zero only if 7 > 7 c, so 


there is a threshold stress for the dislocation nucleation (see (Berdichevsky 


and Le, 2007)). Besides, the existence of the dislocation-free zone in the 


middle of the cross-section can be established. 


7 = 0.01 



Figure 4: The plastic slip /3(r): i) 7 = 0.001, ii) 7 = 0.005, iii) 7 = 0.01. 


For an arbitrary cross section and for ki 7 ^ ^2 the problem does not admit 
exact analytical solution. However, based on the character of solution that 
changes strongly only in the normal direction to the boundary observed in 
the previous case, we may use the asymptotic method to find the solution 
in the thin boundary layer. Take for example the rectangular cross-section 
{xi,X 2 ) e (0,fF) X (0,if), with W and H being its width and height. In 
this case let us assume that there are two boundary layers near the left and 


23 














right boundaries Xi = Q and Xi = hh, and two other boundary layers near 
the bottom and top boundaries X 2 = 0 and X 2 = H. Near the boundaries 
parallel to the X 2 axis the derivative with respect to X 2 can be neglect as 
compared to the derivative with respect to Xi, while the opposite is true 
near the boundaries parallel to xi-axis. In the middle of the cross section 
the plastic slip remains constant. Functional (47) reduces to the sum of four 
integrals, and all of them do not depend on 1. Consider for instance the 
integral along the left boundary layer 


h[P{xi)] = [^(7 - + ]^kil3l]dxi, 

with A being still an unknown length. We omit here the integration over X 2 
(because in this boundary layer it plays just the role of a parameter) and try 
to hnd the minimum among (3. The standard variational calculus leads to 
the Euler equation 

7 - /3 + /ci/3,11 = 0 (49) 

which must be subjected to the boundary conditions 


m = 0 . 


(50) 


The solution of (49) and (50) that does not grow exponentially as xi —)■ cxd 
reads 

/?(xi) = 7(l-e-"^/^). 


Due to the mirror symmetry of the problem, the plastic slip in the boundary 
layer near xi = W must be 


/3(x) =y(l-e-('^-^i)/^). 


Similarly, the solution in the boundary layers parallel to the Xi axis equals 

I 7(1 — e '' 

Thus, the width of the boundary layers parallel to the Xi axis must be of 
the order while that of the boundary layers parallel to the a: 2 -axis must 
be of the order y/ki. Since I depends only on the normal coordinate to 
the boundary, the dislocation lines must be parallel to the boundary of the 
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cross section except at four corners where the asymptotic solution becomes no 
longer valid (the numerical solution should lead to a smoothing of the corners 
of dislocation loops). We see that near the vertical boundaries the dislocation 
loops have the edge character, while near the horizontal boundaries they 
have the screw character. This agrees well with the widths of the boundary 
layers determined by the corresponding contributions of the edge and screw 
components to the energy of the dislocation network. 

For very thin rectangular cross section with H <^W we may neglect the 
influence of the edges near xi = 0 and xi = Why considering the dislocation 
network in the central part of the beam. In this case we will have only screw 
dislocations which pile up against two obstacle at X 2 = 0 and X 2 = H. The 
solution (51) reduces to that found in (Berdichevsky and Le, 2007). 


7. Conclusion 


In this paper we have developed the nonlinear CDT for crystals contain¬ 
ing curved dislocations based on the LEDS-hypothesis. The completely new 
set of equilibrium equations, boundary conditions and constitutive equations 
have been derived from the principle of minimum free energy. As the out¬ 
come, we have obtained the system of strongly nonlinear partial differential 
equations for the placement, the plastic slip, and the loop function. In the 
case of non-vanishing resistance to dislocation motion we have derived the 
governing equations from the variational equation that takes the dissipation 
into account. We have extended the theory to the case of multiple slip and 
simplihed it for small strains. The application of the theory has been il¬ 
lustrated on the problem of single crystal beam having one slip system and 
deforming in simple shear. Under the simplihed assumption ki = /c 2 the ana¬ 
lytical solution of this problem has been found for the circular cross section. 
For arbitrary cross sections the problem has been solved by the asymptotic 
method. We have shown that the asymptotic solution found for the rectan¬ 


gular cross section reduces to the well-known solution in (Berdichevsky and 


Le, 2007) if it is thin and long. 


Acknowledgments 

The hnancial support by the German Science Foundation (DFG) through 
the research projects LE 1216/4-2 and GPOl-G within the Gollaborative 
Research Genter 692 (SFB692) is gratefully acknowledged. 


25 










References 


Acharya, A., 2001. A model of crystal plasticity based on the theory of con¬ 
tinuously distributed dislocations. J. Mech. Phys. Solids 49, 761-784. 

Acharya, A., Bassani, J.L., 2000. Lattice incompatibility and a gradient the¬ 
ory of crystal plasticity. J. Mech. Phys. Solids 48, 1565-1595. 

Berdichevsky, V.L., 2006a. Continuum theory of dislocations revisited. Con¬ 
tinuum Mech. Therm. 18, 195-222. 

Berdichevsky, V.L., 2006b. On thermodynamics of crystal plasticity. Scripta 
Materialia 54, 711-716. 

Berdichevsky, V.L., Le, K.C., 2007. Dislocation nucleation and work hard¬ 
ening in anti-planed constrained shear. Continuum Mech. Therm. 18, 455- 
467. 

Berdichevsky, V.L., Sedov, L.I., 1967. Dynamic theory of continuously dis¬ 
tributed dislocations. Its relation to plasticity theory. Appl. Math. Mech. 
(PMM) 31, 989-1006. 

Bilby, B.A., Gardner, L.R.T., Stroh, A.N., 1957. Continuous distributions 
of dislocations and the theory of plasticity, in: Extrait des actes du JX® 
congres international de mecanique appliquee, pp. 35-44. 

Engels, P., Ma, A., Hartmaier, A., 2012. Continuum simulation of the evolu¬ 
tion of dislocation densities during nanoindentation. Int. J. Plasticity 38, 
159-169. 

Curtin, M.E., 1981. An Introduction to Continuum Mechanics. Academic 
Press, New York. 

Curtin, M. E., 2002. A gradient theory of single-crystal viscoplasticity that 
accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 
5-32. 

Curtin, M.E., Anand, L., Lele, S.P., 2007. Gradient single-crystal plasticity 
with free energy dependent on dislocation densities. J. Mech. Phys. Solids 
55, 1853-1878. 


26 



Hansen, N., Knhlmann-Wilsdorf, D., 1986. Low energy dislocation strnctnres 
dne to nnidirectional deformation at low temperatnres. Mater. Sci. Eng. 
81, 141-161. 

Hochrainer, T., Sandfeld, S., Zaiser, M., Gnmbsch, P., 2014. Continnnm 
dislocation dynamics: towards a physical theory of crystal plasticity. J. 
Mech. Phys. Solids 63, 167-178. 

Hochrainer, T., Zaiser, M., Gnmbsch, P., 2007. A three-dimensional contin¬ 
nnm theory of dislocation systems: kinematics and mean-held formnlation. 
Phil. Mag. 87, 1261-1282. 

Hnghes, D.A., Hansen, N., 1997. High angle bonndaries formed by grain 
snbdivision mechanisms. Acta Mater. 45, 3871-3886. 

Kaluza, M., Le, K.G., 2011. On torsion of a single crystal rod. Int. J. Plas¬ 
ticity 27, 460-469. 

Kochmann, D.M., Le, K.G., 2008a. Dislocation pile-nps in bicrystals within 
continnnm dislocation theory. Int. J. Plasticity 24, 2125-2147. 

Kochmann, D.M., Le, K.G., 2008b. Plastic deformation of bicrystals within 
continnnm dislocation theory. Math. Mech. Solids 14, 540-563. 

Kochmann, D.M., Le, K.G., 2009. A continnnm model for initiation and 
evolntion of deformation twinning. J. Mech. Phys. Solids 57, 987-1002. 

Koster, M., Le, K.G., Ngnyen, B.D., 2015. Formation of grain bonndaries in 
dnctile single crystals at hnite plastic deformations. Int. J. Plasticity 69, 
134-151. 

Kroner, E., 1992. Mikrostrnktnrmechanik. GAMM-Mitteilnngen 15, 104-119. 

Knhlmann-Wilsdorf, D., 1989. Theory of plastic deformation: - properties of 
low energy dislocation strnctnres. Mater. Sci. Eng. A 113, 1-41. 

Knhlmann-Wilsdorf, D., 2001. Q: Dislocations strnctnres - how far from eqni- 
librinm? A: Very close indeed. Mater. Sci. Eng. A 315, 211-216. 

Knhlmann-Wilsdorf, D., Hansen, N., 1991. Geometrically necessary, inciden¬ 
tal and snbgrain bonndaries. Scripta Metall. Mater. 25, 1557-1562. 


27 



Laird, C., Charsley, P., Mughrabi, H., 1986. Low energy dislocation struc¬ 
tures produced by cyclic deformation. Mater. Sci. Eng. 81, 433-450. 

Le, K.C., 2010. Introduction to Micromechanics. Nova Science, New York. 

Le, K.C., Giinther, C., 2014. Nonlinear continuum dislocation theory revis¬ 
ited. Int. J. Plasticity 53, 164-178. 

Le, K.C., Nguyen, B.D., 2012. Polygonization; Theory and comparison with 
experiments. Int. J. Eng. Sci. 59, 211-218. 

Le, K.C., Nguyen, B.D., 2013. On bending of single crystal beam with con¬ 
tinuously distributed dislocations. Int. J. Plasticity 48, 152-167. 

Le, K.C., Sembiring, P., 2008a. Analytical solution of plane constrained shear 
problem for single crystals within continuum dislocation theory. Arch. 
Appl. Mech. 78, 587-597. 

Le, K.C., Sembiring, P., 2008b. Plane-constrained shear of a single crystal 
strip with two active slip-systems. J. Mech. Phys. Solids 56, 2541-2554. 

Le, K.C., Sembiring, P., 2009. Plane constrained uniaxial extension of a single 
crystal strip. Int. J. Plasticity 25, 1950-1969. 

Mayeur, J.R., McDowell, D.L., 2014. A comparison of Gurtin type and mi- 
cropolar theories of generalized single crystal plasticity. Int. J. Plasticity 
57, 29-51. 

Oztop, M.S., Niordson, G.F., Kysar, J.W., 2013. Length-scale effect due to 
periodic variation of geometrically necessary dislocation densities. Int. J. 
Plasticity 41, 189-201. 

Ortiz, M., Repetto, E.A., 1999. Nonconvex energy minimization and disloca¬ 
tion structures in ductile single crystals. J. Mech. Phys. Solids 47, 397-462. 

Ortiz, M., Repetto, E.A., Stainier, L., 2000. A theory of subgrain dislocation 
structures. J. Mech. Phys. Solids 48, 2077-2114. 

Sandfeld, S., Hochrainer, T., Zaiser, M., Gumbsch, P., 2011. Gontinuum 
modeling of dislocation plasticity: Theory, numerical implementation, and 
validation by discrete dislocation simulations. Journal of Materials Re¬ 
search 26, 623-632. 



Sandfeld, S., Thawinan, E., Wieners, C., 2015. A link between microstructure 
evolution and macroscopic response in elasto-plasticity: formulation and 
numerical approximation of the higher-dimensional continuum dislocation 
dynamics theory. Int. J. Plasticity 72, 1-20. 

Sedov, L.I., 1965. Mathematical methods of constructing models of contin¬ 
uum media. Usp. Matem. Nauk. 20, 123-182. 

Wulhnghoff, S., Bohlke, T., 2015. Gradient crystal plasticity including 
dislocation-based work-hardening and dislocation transport. Int. J. Plas¬ 
ticity 69, 152-169. 

Zhu, X., Dai, S., Xiang, Y., 2013. Numerical simulation of dynamics of dislo¬ 
cation arrays and long-range stress helds of nonplanar dislocation arrays. 
Int. J. Plasticity 43, 85-100. 

Zhu, X., Xiang, Y., 2014. Continuum framework for dislocation structure, 
energy and dynamics of dislocation arrays and low angle grain boundaries. 
J. Mech. Phys. Solids 69, 175-194. 


29 



